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ABSTRACT 

Using data obtained by the DIRBE instrument on the COBE spacecraft, we 
present the mean 3.5 — 240 fim spectrum of high latitude dust. Combined with 
a spectrum obtained by the FIRAS instrument, these data represent the most 
comprehensive wavelength coverage of dust in the diffuse interstellar medium, spanning 
the 3.5 — 1000 /xm wavelength regime. At wavelengths shorter than ~ 60 the 
spectrum shows an excess of emission over that expected from dust heated by the local 
interstellar radiation field and radiating at an equilibrium temperature. The DIRBE 
data thus extend the observations of this excess, first detected by the IRAS satellite at 
25 and 12 ^m, to shorter wavelengths. The excess emission arises from very small dust 
particles undergoing temperature fluctuations. However, the 3.5-to-4.9 /im intensity 
ratio cannot be reproduced by very small silicate or graphite grains. The DIRBE data 
strongly suggest that the 3.5 - 12 fim emission is produced by carriers of the ubiquitous 
3.3, 6.2, 7.7, 8.6, and 11.3 /mi solid state emission features that have been detected in 
a wide variety of astrophysical objects. The cariers of these features have been widely 
identified with polycyclic aromatic hydrocarbons (PAHs). 

Our dust model consists of a mixture of PAH molecules and bare astronomical 
silicate and graphite grains with optical properties given by Draine &; Lee. We obtain 
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a very good fit to the DIRBE spectrum, deriving the size distribution, abundances 
relative to the total hydrogen column density, and relative contribution of each dust 
component to the observed IR emission. At wavelengths above 140 fim the model is 
dominated by emission from T ~ 17 — 20 K graphite, and 15 — 18 K silicate grains. 
The model provides a good fit to the FIRAS spectrum in the 140 — 500 /xm wavelength 
regime, but leaves an excess Galactic emission component at 500 — 1000 |im. The 
nature of this component is still unresolved. 

We find that (C/H) is equal to (7.3±2.2) 10~ 5 for PAHs, and equal to (2.5±0.8) 10" 4 
for graphite grains, requiring about 20% of the cosmic abundance of carbon to be 
locked up in PAHs, and about 70% in graphite grains (we adopt (C/H)q = 3.6 10~ 4 ). 
The model also requires all of the available magnesium, silicon, and iron to be locked 
up in silicates. The power emitted by PAHs is 1.6 10~ 31 W/H atom, by graphite 
grains 3.0 10 -31 W/H atom, and by silicates 1.4 10 -31 W/H atom, adding up to a 
total infrared intensity of 6.0 10 -31 W/H atom, or ~ 2 L & /M & . 

The [C II] 158 //m line emission detected by the FIRAS provides important 
information on the gas phase abundance of carbon in the diffuse ISM. The 158 fj,m line 
arises predominantly from the cold neutral medium (CNM), and shows that for typical 
CNM densities and temperatures C+/H = (0.5 - 1.0) 10~ 4 , which is ~ 14 - 28% of 
the cosmic carbon abundance. The remaining carbon abundance in the CNM, which 
must be locked up in dust, is about equal to that required to provide the observed IR 
emission, consistent with notion that most (;> 75%) of this emission arises from the 
neutral component of the diffuse ISM. 

The model provides a good fit to the general interstellar extinction curve. However, 
at UV wavelengths it predicts a larger extinction. The excess extinction may be the 
result of the UV properties adopted for the PAHs. If real, the excess UV extinction 
may be accounted for by changes in the relative abundances of PAHs and carriers of 
the 2200 A extinction bump. 



Subject headings: interstellangrains, infrared: radiation - spectrum, ISM: Dust, 
Extinction 



1. INTRODUCTION 

A surprising result of the Infrared Astronomical Satellite (IRAS) all-sky survey was the 
detection of a "new" high-latitude diffuse infrared emission component dubbed "cirrus" because 
its morphology resembled that of filamentary clouds in our own atmosphere ( [Low et al. 1984 ). 



This IR emission component was not entirely new, and had been previously detected through 



scattered light at optical wavelengths since the 1970's (e.g. le Vaucouleurs & Freeman 1972, 



Sandage 1976). The surprising aspect of the infrared (IR) observations was their spectra which 
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displayed an excess of emission at 12 and 25 ^m over that expected from dust particles radiating 
at the equilibrium temperature and heated by the local interstellar radiation field (LISRF). Draine 
&; Anderson (1985) attributed this excess to the effect of temperature fluctuations in small grains. 
The grain size distribution given by Mathis, Rumple, &: Nordsieck (1977; hereafter MRN) was 
therefore required to be extended to very small grain sizes of ~ 3 A ( |Draine &: Anderson 198*5 , 



Weiland et al. 1986 ). The cirrus spectrum may, however, arise from the same population of dust 
particles that give rise to the "unidentified" IR emission features at 3.3, 6.2, 7.7, 8.6, and 11.3 /im. 
These features were detected in a wide variety of astrophysical objects, reflection nebulae, H II 



regions, and planetary nebulae (see Cohen et al. 198"6| and references therein), and the large-scale 



emission in the 3.3 /im band was detected by Giard et al. (1988, 1994). However, so far there 
has been no direct detection of these features in high latitude cirrus clouds, a situation that may 
change when analysis of data from the ISO and IRTS satellite missions are completed. 

The carriers of these features are most commonly identified with a class of molecules called 



polycyclic aromatic hydrocarbons, or PAHs (Duley & Williams 1981, Leger & Puget 1984 



Allamandola et al. 1985, Puget et al. 1985| ) . However, that identification is not unique, and 



various other compounds such as hydrogenated amorphous carbon (HACs; e.g. Duley 1989, and 
references therein), or quenched carbonaceous composites (QCCs; e.g. Sakata & Wada 1989, 
and references therein) have been suggested as alternative carriers of these features. Specifically, 
Papoular et al. (1989) suggested demineralized coal or vitrinite for the carriers. The optical 
constants of 10 — 100 A coal grains are such that the grains can maintain equilibrium temperatures 
with near-IR colors in the 500 — 1000 K range, provided that they are exposed to radiative fluxes 
of the order of 10 -2 — 1 W m~ 2 ( puillois et al. 1994 ). These radiative fluxes are a few orders 



of magnitude larger than those of the LISRF, suggesting that coals are not likely carriers of the 
emission features in the diffuse interstellar medium (ISM). Throughout our analysis we will assume 
that PAHs are the carriers of the unidentified IR emission bands. 

Observations of cirrus made by the Diffuse Infrared Background Experiment (DIRBE) and 
Far Infrared Absolute Spectrophotometer (FIRAS) experiments on board the Cosmic Background 
Explorer ( COBE) satellite (|Boggess et al. 1992[) represent the first spectra of these objects 
spanning the entire 3.5 to 1000 /im wavelength region. They therefore provide the strongest 
constraints on the nature of the dust in these regions. Bernard et al. (1994, 1996) presented 
the DIRBE observed cirrus spectrum for low Galactic latitude {\b\ < 10°). In this paper we 
complement their studies by deriving an average high-latitude cirrus spectrum and extending the 
wavelength coverage to 1000 fj,m. We also present the spectrum of select individual high-latitude 
cirrus clouds in order to examine the spectral variations in these objects. 

In general, the determination of the cirrus spectrum from the DIRBE maps is complicated as 
a result of the presence of strong foreground emissions from the interplanetary dust cloud and at 
3.5 and 4.9 yum from Galactic stellar emission. In §2 we describe the data set used in the analysis, 
and briefly reiterate the method developed by Arendt et al. (1996; also presented by [Weiland^ 



Arendt, &; Sodroski 1996 ) for the extraction of the cirrus spectrum. The IRAS data provided clear 
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evidence for the presence of transiently heated dust particles in cirrus clouds flDraine fc Anderson 
1985D - However, as a result of the limited wavelength coverage of the IRAS (A = 12, 25, 60, and 
100 /im), the cirrus data could not constrain the composition of the very small dust particles. In 
§3 we show that any extension of the MRN size distribution to very small grains is inconsistent 
with the DIRBE data. Adopting the PAH model described by Desert, Boulanger, & Puget (1990) 
to characterize the carriers of the 3.3 — 11.3 /mi features (see Appendix A. 2), we show that the 
DIRBE 3.5 to 25 /mi data are consistent with the predicted spectrum from these macro- molecules. 
A detailed description of the dust model used in this paper is presented in §4, and in Appendix 
A.l we describe the method used to calculate the PAHs heat capacities used in this paper. In §5 
we describe the fitting procedure, and discuss the uniqueness of the various model parameters 
derived from fitting the DIRBE observations. We also calculate the abundances of the various dust 
constituents relative to hydrogen, and calculate the extinction predicted by the model, comparing 
it to the observed average interstellar extinction curve. The [C II] 158 /im line emission detected 
by the FIRAS instrument is used to calculate the gas phase abundance of carbon in the diffuse 
ISM (Appendix A. 3), which provides an independent estimate of the amount of carbon required to 
be locked up in dust. In §6 we compare the far-IR spectrum predicted by the model to the FIRAS 
observations. The results of our paper are briefly summarized in §7. 



2. THE DERIVATION OF THE CIRRUS SPECTRUM 

A description of the derivation of the dust spectrum in the general ISM was presented by 
Weiland, Arendt, & Sodroski (1996), based on the detailed work of Arendt et al. (1996). The data 
used in the derivation consists of the CO-RE/DIRBE sky maps of emission at A = 1.25, 2.2, 3.5, 
4.9, 12, 25, 60, 100, 140, and 240 /mi (with corresponding bandwidths of 570, 225, 220, 82.1, 135, 
41.0, 23.2, 9.74, 6.17, and 4.96, in units of 10 11 Hz), from which the contributions of interplanetary 
dust particles and of discrete bright stellar and extragalactic sources and unresolved Galactic 



starlight has been removed (Kelsall et al. 1996, Arendt et al. 1996). The remaining intensity 



consists therefore of IR emission from Galactic interstellar dust, and from a diffuse extragalactic 
background which is assumed to be spatially isotropic. 

The spectrum of the ISM, relative to the 1(100 /im) intensity, was obtained by deriving the 
slope of the 1(A) versus 1(100 /im) correlation for the high latitude regions of the sky that are 
unaffected by residual interplanetary dust emission. The slope of the correlation is insensitive 
to any isotropic terms, and the spectrum is therefore free of any extragalactic or any isotropic 
Galactic emission components. The spectrum was then scaled to the Galactic H I column density 
by correlating the 7(100 /im) intensity to N(H I) over the same regions of the sky. However, by 
construction, the resulting IR emssion does not originate only from the H I gas, but from all gas 
phases that correlate with the 1(100 /im) intensity. 

At wavelengths A < 4.9 /im the subtraction of the stellar emission component is not sufficiently 
accurate to bring out the dust emission component in the residual intensity maps. Therefore, 
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Arendt et al. used near-IR reddening-free colors to identify this component using DIRBE 1.25, 
2.2, 3.5, and 4.9 /jm data. The procedure is equivalent to identifying dust emission components in 
color-color plots as the points that are offset from the stellar reddening curve. The method could 
only be applied to Galactic latitudes with \ b\ < 25° because of the presence of residual artifacts 
from the interplanetary dust emission at higher latitudes. As a result there may be a systematic 
offset between the the 3.5 and 4.9 /mi intensities and the rest of the spectrum, if the ratio of the 
3.5 or 4.9 to 100 /mi intensity is a function of latitude. 

The cirrus spectrum at wavelengths longward of 100 /mi was derived from the FIRAS Pass 3 
data (see Fixsen 1994) from which the microwave background and the dipole have been subtracted, 
in a fashion similar to that used with the DIRBE data at A > 12 /mi. The first step in this 
derivation consisted of degrading the Galactic component of the DIRBE 100 /mi emission to the 
FIRAS 7° beam resolution using the actual FIRAS beam pattern derived from observations of 
the moon. The full FIRAS high frequency data set was used for 100 < A(/xm) < 500, while the 
combined low frequency data set was used for the 500 — 5000 /mi wavelength region. The spatial 
correlation with the DIRBE 100 /mi intensity map was performed for each of the 167 FIRAS 
frequency channels. In both frequency regimes, the data were limited to |6| > 45° and weighted by 
the FIRAS weights. In the resulting spectrum, the FIRAS intensity matches that of the DIRBE 
240, 140, and 100 /mi in the respective wavelength bands, confirming the consistency of the data 
sets and the derivation procedure. Only the 140 — 1000 /xm spectrum of the FIRAS is presented 
here since the data are too noisy in the remaining frequency channels. 

The derived 3.5 — 1000 /mi cirrus spectrum is therefore an average spectrum of high-latitude 
interstellar dust. It represents the average emission from dust integrated over all the non- isotropic 
ISM components. Dust residing in H II as well as H2 gas can therefore contribute to the emission. 
To examine the gas phase dependence of the dust emission we derived the IR spectra of two 
distinct high latitude clouds. The first cloud (referred herafter as Cloud 1) corresponds to cloud 



numbers MBM 53, 54, 55 in the catalog of high-latitude molecular cirrus clouds (Magnani, Blitz, 



& Mundy 1985). It represents a cloud in which the IR emission is dominated by dust residing in 
the molecular gas and its IR emission is modeled in detail by Dwek, Arendt, Fixsen, & Reach 
(1996). The second (Cloud 2) was referred to as Cloud A by Low et al. (1984) and modeled by 
Draine &; Anderson (1985). This cloud is notable for having the highest column density of any 
feature in the b < —60° H I map of Heiles (1975). 

Figure 1 depicts the various cirrus spectra derived in this paper. The diamonds represent the 
average ISM spectrum, and the spectra of Clouds 1 and 2 are represented by triangles and stars, 
respectively. The average ISM spectrum is normalized to a value of 0.7 MJy sr~ l /iVjy(10 20 cm 2 ) 
at A = 100 /mi. The other spectra are offset by an arbitrary factor for sake of clarity. Uncertainties 
in the data are determined by the RMS variation of the I(A)//(100 /mi) colors in different regions 
of the sky, and equal to 20% (1 d). We note that the average ISM spectrum presented in the 
figure is very similar to that derived by Bernard et al. (1994, 1996). 
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Infrared lines can contribute to the observed emission in the broad DIRBE bands. The 
most important potential contributor in the diffuse ISM is line emission from [C II] 158 /xm in 
the 140 /im DIRBE band. Our calculations (see §5.2 below) show that this line contributes at 
most ~ 2% to the emission in this band. We therefore attribute the observed spectra to thermal 
emission from interstellar dust. 

Table 1 summarizes the various spectra derived in this paper. Also presented in the table 
are the dust color temperatures determined from the flux ratios in adjacent bands. Refering to 
the average cirrus spectrum, the 140-to-240 /mi and the 60-to-100 /mi flux ratios correspond to 
color temperatures of ~ 19 K, and ~ 23 K, respectively (for a A~ 2 emissivity law), whereas the 
12-to-25 (im and the 3.5-to-4.9 /jm flux ratios give color temperatures of ~ 210 K, and ~ 550 K, 
respectively. 

The excess emission at shorter wavelengths has been attributed to the fact that the cirrus 
emission arises from two distinct populations of dust particles (Draine fc Anderson 1985Q . The first 



consists of particles radiating at a relatively narrow range of equilibrium temperatures between 
16 and 20 K, and the second consists of dust particles small enough to undergo temperature 
fluctuations. Their contribution to the emission is first noticable at 60 /mi, as manifested in the 
constancy of the 60-to-100 /mi Galactic IRAS flux ratio with galactic longitudes ( podroski et al 



1989| , podroski et al. 1994| ), and increases at shorter wavelengths to dominate the emission at 
wavelengths below 25 /mi. Comparison of the spectra depicted in Figure 1, or listed in Table 1, 
shows variations in the long and short wavelength color temperatures. These variations probably 
reflect differences in the particle size distributions and/or variations in the intensity or spectrum 
of the interstellar radiation field. An excellent review of the conditions leading to the stochastic 
heating of interstellar dust particles and the observational effects of the resulting temperature 
fluctuations is given by Aannestad (1989). 



3. THE DIRBE EVIDENCE FOR THE PRESENCE OF PAHs 

The IRAS observations of the diffuse ISM were of limited wavelength coverage and could 
therefore not constrain the composition of the very small dust particles. Models by Draine & 
Anderson (1985) and Weiland et al. (1986) could fit the IRAS fluxes by simply extending the size 
distribution in the interstellar graphite-silicate MRN dust model to very small grain sizes (a ~ 3 
A ). Subsequent calculations ( puhathakurta fc Draine 1989 ) have shown that graphite and silicate 



grains with radii less than ~ 4 A, and ~ 5 A, respectively evaporate after residing for ~ 4 x 10 5 
yrs in the LISRF. The extension to grain sizes below these values is therefore unrealistic in the 
absense of any replenishing mechanism. Similarly, only PAHs with more than ~ 20 carbon atoms 
will survive evaporation in the LISRF ( pmont 1986| ). 

The DIRBE has more extensive wavelength coverage than the IRAS, and any extension of 
the standard MRN bare silicate-graphite model to small grain sizes that survive evaporation by 
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the LISRF is inconsistent with the DIRBE data. However, the DIRBE 3.5 to 25 data are 
consistent with the predicted spectrum from PAHs, if we adopt the model described by Desert, 
Boulanger, & Puget (1990) to characterize the carriers of the 3.3 — 11.3 \xm features (see Appendix 
A. 2). The "detection" of PAHs with the DIRBE is demonstrated in Figures 2a— 2d which compare 
various observed intensity ratios with calculated in-band DIRBE intensity ratios of individual 
grain sizes for various grain compositions. Figure 2a shows this comparison for the 3.5-to-4.9 fim 
intensity ratio (hereafter designated as #3.5/4.9). The figure illustrates that no size distribution 
of silicate or graphite grains can reproduce the observed intensity ratio, since even the smallest 
grains (a = 4 A, for graphite, and a = 5 A for silicates) have intensity ratios that fall below the 
observed one. The value of #3.5/4.9 is expected to decrease with increasing grain size, since the 
amount of short wavelength emission depends strongly the maximum temperature reached by the 
stochastically-heated particles, and for a given radiation field, this temperature decreases with 
increasing grain size. On the other hand, small PAHs can produce larger values of #3. 5/4.9 because 
they posses a strong emission feature at 3.3 fj,m. The figure also gives a good indication of the 
range of PAH sizes that is required to fit the observed intensity ratio. PAHs are two-dimensional 
structures, and the relation between the PAH radius and its number of carbon atoms, N c , is 
given by a(A)=0.913\/^Vc (Desert, Boulanger, & Puget 1990). The figure suggests an average 
PAH radius of ~ 7 A, which consists of about 60 carbon atoms. This average PAH size is in 
reasonable agreement with the value of ~ 90 carbon atoms estimated by Leger, d'Hendecourt, & 
Defourneau (1989) to comprise a typical PAH molecule giving rise to the observed 3.3 to 11.3 ^m 
color temperature. Figure 2b illustrates the same effect for #12/25) the I (12 fim)/ 1 (25 fj,m) intensity 
ratio. In this case, very small graphite and silicate grains can reproduce the observed intensity 
ratio, illustrating why a simple extension of the standard MRN dust model to very small grain 
sizes could explain the IRAS data. Figures 2c and 2d continue this comparison for #60/100 an d 
#140/240) respectively, both covering wavelength regimes where dust radiating at the equilibrium 
temperature dominates the emission. This fact is well illustrated in these figures which show that 
PAHs have significantly higher color temperatures than those implied by the observations. PAHs 
therefore can not be important contributors to the emission at these wavelengths. 



4. A DESCRIPTION OF THE DUST MODEL 

The IR flux from any given line of sight is the sum of emissions from dust residing in all the 
diffuse ISM phases along that line of sight. The specific IR intensity from a population of dust 
particles characterized by a size distribution f(a) in the radius interval a m in < a < a ma x, and a 
dust-to-gas mass ratio Z^, can therefore be written in the most general form as: 

/(A) = x N H I"™* da f(a) a(a, A) / dT V(a, T) B X (T) (1) 

<m d > J amin J 

where \i is the mean atomic weight of the gas in amu, Nh = Nh i + Nh ii + 2iV# 2 is the total 
H-column density along the line of sight, < m<± >= !aZZ ^ a f( a )lfP a3 i * s * ne size-averaged mass 
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of the dust population, p is the mass density of a dust particle, a(a, A) is its cross section at 
radius a and wavelength A, B\(T) is the Planck function, and V(a,T)dT is the probability that 
the temperature of a dust particle will be between T and T + dT. For dust particles radiating at 
the equilibrium dust temperature, T eq , the function V(a,T) is simply a 5— function at T = T eq . 
Implicit in equation (1) is a summation over the various grain compositions. A similar equation 
applies for PAHs, except that we represent the PAH size distribution as a function of the number 
of carbon atoms, rather than grain radius. 

The dust model used in this paper consists of a mixture of bare silicate and graphite particles, 
and PAH molecules. An important physical property of the radiating dust particles is their heat 
capacity. This quantity determines the temperature fluctuations, and hence the function V(a,T) 
of the various dust particles. Appendix A.l describes the heat capacities used for the various grain 
constituents, with emphasis on PAHs, for which we used the group additive method to determine 
their value. The relation between V{a,T) and the various grain parameters and the LISRF 
has been described by Aanestad (1989). The LISRF used in this paper was taken from Mathis, 
Mezger, & Panagia (1983). Also summarized in the Appendix are the radiative cross sections used 
in the model. Graphite and silicate optical properties were calculated from Mie theory using the 
optical constants of Draine & Lee (1984). Optical properties of the PAH were adopted from the 
work of Leger, d'Hendecourt, &; Defourneau (1989), as summarized by Desert, Boulanger, &: Puget 
(1990). 

Given the dust composition and physical and optical properties, a dust model must specify 
the size distribution and relative abundance of its various constituents. The distributions in grain 
radii, a, adopted for the various grain material are: 

(a) for silicates: 

f(a) = a 1 a min <a< a max (2) 

with {a m i n , 7, a max } = {0.0050 fj,m, —3.5, 0.25 /mi}, their MRN value, except that the lower size 
limit was extended from 0.025 /im, the nominal MRN value, to a lower value of 0.0050 /im; 

(b) for graphite: 

f(a) = a 11 a min < a < a b 

= a 12 ab <a< a max (3) 

allowing for a break in the slope of the power law at some intermediate radius aj,; 
and (c) for PAHs: 

The grain size distribution is written as a power law in N c , the number of carbon atoms in a PAH 
molecule: 

f(N c ) = N? N cl <N C < N c2 (4) 

The number of carbon atoms in the smallest PAH was fixed at N c \ = 20 since smaller PAHs will 
evaporate when exposed to the LISRF ( pmont 1986 ). The number of carbon atoms characterizing 
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the largest PAH, and the power law of the PAH size distribution, 7 P , were taken to be free 
parameters of the model. 

We have considered two models using with different constraints on the graphite size 
distribution. For both models the graphite size distribution is characterized by {72, a max } = 
{—3.5, 0.25 /iin} for radii above a&. The models differ in their characterization of the distribution 
of very small graphite grains. The first model (model A) considers the very small graphite 
grains and the PAHs to be independent populations of particles. We adopt the MRN grain size 
distribution for this model, but allow a m i n , the lower limit of the size distribution, to attain values 
below the nominal MRN value of 0.0050 /im. In this model, 71 = 72 = —3.5, the MRN value. 

In the second model (model B), the graphite size distribution is smoothly joined to that of the 
PAHs, a distribution one would expect if PAHs were created from the destruction and processing 
of graphite grains in the ISM. The size of the smallest graphite grain a m i n was therefore chosen so 
that the number of carbon atoms in a spherical graphite particle of radius drain 

would equal 

the number of carbon atoms in the largest PAH, that is, 

aminipm) = 1.29 x lO" 4 N% 3 (5) 

The power law of the size distribution in the {a m i n , a^} size interval was chosen so that dn(N c )/dN c 
is continuous across trie PAH-graphite boundary at Qrnin- 

The value of 71 is therefore related to 

1p by 

71 = 3 7p + 2 (6) 

The constraints on a m i n (graphite) and and on 71 and 7 P as expressed in equations (5) and 
(6) are quite arbitrary, but they have the aesthetic virtue of providing a smooth transition from 
graphite grains to PAHs. 

The upper limit of the size distribution of the graphite and silicate grains is uncertain. 
It cannot be constrained from extinction measurements, since large grains are essentially gray 
particles. However, the far-IR emission is sensitive to the value of the upper limit, since larger dust 
particles are colder. Unfortunately, it is difficult to disentangle the effects of grain size from the 
effects of the intensity of the interstellar radiation field. A higher value of the LISRF will increase 
the average dust temperature, an effect that can also be achieved by reducing the value of a max 
by an appropriate value. With this ambiguity in mind, we simply followed previous dust models 



(e.g. Draine fc Anderson 1985 ), and chose the value of this cutoff to be at a max = 0.25 /j,m. 



To summarize, the dust model is characterized by 14 parameters: 4 for silicates 
Z si i}; 6 for graphite {a m j„, 71, a&, 72, a max , Z gr f}; and 4 for the PAHs 
{iVci, 7 p , iV C 2, Zpah}- Of these parameters, 8 are fixed, and 6 are allowed to vary in the 
calculations. Table 2 summarizes the various parameters of the dust model. 



5. MODEL RESULTS 
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5.1. General 

To find the dust model that best fits the observed spectrum we mapped out the 3-dimensional 
X 2 space spanned by the model parameters: {7^, N C 2, a m in (graphite)}, or {7^, N C 2, at>} depending 
on the characterization of the graphite grain size distribution (model A or B, respectively). For 
each value of these 3 parameters, the PAH-to-gas mass ratio, Zpah was determined from fitting 
DIRBE 3.5, 4.9, and 12 /im observations with the PAH spectrum. The graphite and silicate 
dust-to-gas mass ratios, Z gr f and Z s u, respectively, were then determined from a x 2_ fit of the 
model to the residual spectrum. 

We expect some correlation between the various dust parameters. For example, N C 2 and 
7 p should be strongly correlated as suggested by Figure 2a. The figure shows that the observed 
3.5-to-4.9 fim flux ratio is about equal to that produced by PAHs with ~ 60 carbon atoms. The 
ratio can however be equally well reproduced by a range of PAH sizes, from N& = 20 to some 
upper limit N C 2- The larger the value of N& the less these PAHs will contribute to the 3.5 p 
emission, and the more the spectrum will need to weighted towards the smaller PAH sizes. Large 
values of N C 2 will therefore correspond to steep power laws, i.e. large values of |7 P |. 

For a given spectrum, grain abundances are primarily determined by their opacities per unit 
mass. These quantities are relatively insensitive to grain radius, so that the derived abundances 
are only sensitive to details of the IR spectrum. For example, the abundance of silicates is sensitive 
to the 140-to-240 /im flux ratio. Exposed to the LISRF, 0.20 /im radius silicate grains attain a 
lower equilibrium temperature (T eq ~ 15 K) than identically sized graphite particles for which 
T eq 18 K. Smaller values of .R140/240; wm require colder dust to participate in the emission, 
and will therefore lead to larger values of Z s n. 

Table 3 lists the values of the parameters of four models for the average ISM that produced the 
lowest values of x 2 - For each choice of the graphite grain size distribution (model A or B in Table 
2), the four models represent the extreme values of N C 2 that still provide a good fit to the data. 
The various models produce almost indistinguishable fluxes with minor differences in the spectra 
of the various grain constituents. Typical values of N C 2 are between ~100 and 170. Figure 3 
presents the spectrum for dust model A with N C 2 = 100. The figure shows that PAHs dominate the 
emission at wavelengths A <J 20/im, whereas stochastically heated graphite grains are important 
contributors to the emission in the narrow wavelength region of 40 <^ A(/im) <^ 70. Graphite and 
silicate grains contribute about equally to the emission at wavelengths above ~ 300 /j,m. The 
main effect of the silicate grains is to broaden the peak of the spectrum around ~ 140 /im. The 
power emitted by PAHs is 1.6 10~ 31 W/H atom. Graphite grains emit 3.0 10 -31 W/H atom, 
and silicates emit 1.4 10~ 31 W/H atom. The total infrared intensity is 6.0 10 -31 W/H atom, or 



1.9 Lq/Mq, compared to an average value of ~ 3.0 Lq/Mq found for the inner Galaxy (Sodroski 



et al. 1994 



PAH size distributions derived here are very similar to that derived by Desert et al. (1990). 
Translating the results from their Table 2 to parameters used in this paper we get for their model: 
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N c i = 19; N C 2 = 170; 7 P = —1.5. Similar results were obtained by Siebenmorgen & Kriigel (1992) 
who derived a typical PAH cluster size of 150 carbon atoms. 



5.2. Dust and PAH Abundances 

Assuming that the H2 and H + contributions to the total H-column density are negligible, we 
get that the total dust-to-gass mass ratio of the model, Z^nst, is equal to 0.0080, consistent with 
the value of 0.0073 allowed from cosmic abundance considerations ( Anders and Grevesse 1989| ). 



The small excess may be attributed to the fact that Nh may be larger than Njj j, which would 
indicate that dust in H I-correlated ionized or molecular medium contributes to the IR emission. 
The size of the discrepancy seems, however, to indicate that this contribution is small compared to 
that of the dust residing in the neutral medium. We also note Z^ ust may be larger than the quoted 
value, since an additional mass of dust may be required to explain the cold dust component in the 
FIRAS spectrum (see §6 for more details), or may be present in the ISM but too cold to manifest 
itself in the current observations. 

Similar results are obtained for the carbon abundance required by the model to be in 
either the PAHs or the graphite grains. The mass fraction of carbon in PAHs is (see Table 3): 
Zpah = 0.00062, (giving C/H = 7.3 10 -5 for a mean molecular weight fi = 1.42), and for graphite 
grains Z gr f = 0.0021, (C/H = 2.5 10 -4 ). The calculated carbon mass fraction in the solid phase 
is thus Zc( so iid) = 0.0027 (C/H = 3.2 10 -4 ), slightly less than the cosmic abundance value of 
Z C (Q) = 0.0030 [ (C/H) Q = 3.6 10~ 4 ] ( [Anders and Grevesse 19891 ) . 

The "missing" carbon must be in the gas phase and singly ionized (see discussion in Appendix 
A. 3). Its abundance can therefore be derived from the [C II] 158 fim line emission detected by 
the FIRAS flWright et al. 1991| , |Bennett et aL 1994j ). The 7(158 fim)/ 1 (100 /im) line ratio is 



latitude dependent, and subtracting the dust continuum emission we find a C + cooling rate per 
H-atom of 1.45 10 -33 W/H atom for the high-latitude ISM. This is about half the average value 
of 2.65 10 -33 W/H atom derived by Bennett et al. (1994) from data that included lower latitude 
observations. The power in the 158 /im line constitutes about 2% of the total energy received in 
the DIRBE 140 /im band, and its intensity implies a [C + ]/[H I] ratio of ~ (0.5 — 1) 10~ 4 (see 
Appendix A. 3), which is about equal to 14 — 28% of the cosmic carbon abundance. The total 
carbon to H-atom ratio is therefore C/H = (3.7 — 4.2) 10~ 4 . This value is consistent with the 
cosmic abundance limit, confirming previous results that any contribution of the ionized gas to 
the IR emission is small (<^ 15%). Assuming that all the IR emission originates from the H I gas, 
we get that ~ 18% of the carbon is locked up in PAHs, ~ 62% in graphite, and the remaining 
~20% is in the gas phase. 

The abundance of PAHs needed to account for the IR emission can also be estimated in the 
following two ways: the first by balancing the energy absorbed by the PAHs from the LISRF to 
the observed emission in the 3.5 — 12 /im wavelength regime, and independently, by calculating the 
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PAH optical depth at some representative wavelength dominated by PAH emission. 

In the first method the fraction of carbon atoms locked up in PAHs is determined by the ratio 
of the energy absorbed per C-atoms locked up in PAHs: 

c J U LIS rf(X) a PAH {\)d\ = 2.6 1(T 27 W/(C atom) (7) 

and the observed IR luminosity emitted by the PAHs per H atom: 

4ttI pah /N h = 1.6 10~ 31 W/(H atom) (8) 

where in eq. (7), Ulisrf is the energy density of the LISRF, and a pah is the PAH cross section 
per carbon atom in the molecule. The resulting value is: 

{C/H) PAH = 6.1 1(T 5 (9) 

similar to the value of (C/H) PAH = 5.4 10" 5 derived by Puget & Leger (1989). 

In the second method we calculate the abundance of PAHs required to account for the 
observed IR emission. These calculations are a simple analytical representation of the detailed 
model used to derive the results presented in Table 3. We first assume that all the intensity at a 
given wavelength A originates from PAHs radiating at an equilibrium temperature T eq . The PAH 
optical depth tp A h per H column density will then be given by: 

tpahW = IW/Nh , . 

N H B x (T eq ) { ' 

where B\(T) is the Planck function. The PAH optical depth can also be expressed as: 

tpah(A) = N c a PAH (X) (11) 

where Nc is the column density of carbon atoms locked up in PAHs. The C/H ratio locked up in 
PAHs is then derived from equations (7) and (8): 

trim iW/Nh , . 

{C/H)PAH ~ a PAH (X) B x (T eq ) (12) 

We now assume that all the emission in the DIRBE 12 ^m is emitted by PAHs. From Table 1 
we get that I±2 ^m/Nn = 3.16 1CP 16 Jy sr~ l cm~ 2 . PAHs contributing to the emission will be at 
temperatures of about 425 K, for which B\2 Mm (425 K) = 1.46 10 15 Jy sr^ 1 , and the PAH cross 
section per unit C-atom is ap A n(l2 fim) = 1.4 1CT 21 cm 2 . The resulting value of (C/H)p A h is 
1.5 1(T 10 . 

A major assumption made in the derivation of this number is that the PAHs radiate at an 
equilibrium temperature, when in fact, they are undergoing temperature fluctuations. Therefore, 
at any given time, the 12 /xm observations sample only a fraction of the PAH population, those 
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that happen to be at temperatures around 425 K. The fraction of "flickering" PAHs is given by 
T cool(T)/ < T a b s >) where t coo i(T) is the cooling time at temperature T, and < r a i, s > is an average 
time between photon absorptions. The latter quantity can be simply calculated from: 

< r abs >^= J {^)U LIS rf{\) a PAH (X)dX (13) 
where h is the Planck constant. The PAH cooling time is given by: 

t m-i 1 4 < d PAH (T) > aT* 

r cool (T) =- ^— (14) 

where < 5pah{T) >= 1.4 10~ 21 cm 2 is the Planck-averaged value of the PAH cross section per C 
atom, and Cpah{T) = 2 W~ w erg K~ l is the PAH heat capacity per C atom given by eq. (A-4). 
The resulting cooling time around T = 425 K is ~ 8 s. The absorption time is size dependent. To 
yield a {C/H) pah value consistent with the value obtained in eq. (9) the fraction of flickering 
PAHs must be equal to 1.5 10~ 10 /6.1 10 -5 = 2.4 10~ 6 , that is, the average time between photon 
absorption should be 3.2 10 6 s. The required PAH size is N c « 80, consistent with the typical 
PAH sizes giving rise to the IR emission. 

The derived fraction of carbon locked up in PAHs is about 2 times the value derived by 
Desert et al. (1990; C/H = 3 10" 5 ), Verstraete & Leger (1992; C/H = 4 10" 5 ), or Siebenmorgen 
& Kriigel (1992; C/H = 3 10~ 5 ), but comparable to that estimated by Puget & Leger (1989; 
C/H = 5.4 10~ 5 ). The abundance of the PAHs is directly proportional to I pah {12 M m )/-^ff) an d 
inversely proportional to the wavelength integrated value of Ulisrf x &pah- Differences between 
the PAH abundance derived here and those derived by Desert et al. can be completely accounted 
for by the fact that the value of the vIpah (12 [xm)/vI{Q0 /im) ratio in our spectrum is larger by 
a factor of ~ 1.45 than their value. Differences between our PAH abundance determination and 
that of Verstraete & Leger are due to their use of larger UV-visual cross-sections for the PAHs. 

Uncertainties in the derived PAH abundance are dominated by uncertainties in the value of 
the LISRF, the ipAff(12 fim)/NH ratio, and the UV-optical cross sections of the PAHs. Total 
uncertainties in the LISRF are about 15% (Mathis, Mezger, & Panagia 1983). The uncertainty 
in the ratio of the 12 /jm intensity to H-column density ratio is determined by the combined 
uncertainties in the 1(100 fimj/Nff and 1(12 /im)/I{W0 fxm) ratios from which it was derived 
(Arendt et al. 1996). Uncertainties in each of these quantities are about 10% and 20%, respectively. 
Combined with the uncertainties in the LISRF, we estimate the uncertainty in the PAH abundance 
to be about 30%. Uncertainties resulting from those in the UV-optical cross section of the PAHs 
are harder to evaluate. In our model we have used a value of ~ 3 10~ 18 cm? /C atom in the 1000 
- 3000 A wavelength range, which is a conservative value according to Allamandola, Tielens, & 
Barker (1989). Significantly larger UV absorption cross sections will reduce the required PAH 
abundance. Laboratory data suggest that PAH absorption cross sections in this wavelength region 
can be larger by about a factor of ten (Joblin, Leger, & Martin 1992). Such large cross sections will 
make them significant contributers to the UV extinction in this wavelength regime, and potential 
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carriers of the 2200 A extinction feature. This possibility will require significant modifications to 
current dust models, and will not be further considered here. In summary, for the UV-optical 
cross sections adopted in this paper, the estimated uncertainty in the PAH abundance is about 
30%. The same uncertainty applies to the silicate and graphite abundances derived in this paper. 

There has recently been an accumulating body of evidence suggesting that the carbon 
abundance in the gas phase of the local ISM is (C/H)jsm = (2.25 ± 0.5) 10~ 4 (Snow & Witt 
1995), about 60% of its solar value. This lower value severely constrains the amount of carbon 
that can be locked up in dust, leading various authors (Kim & Martin 1996, Mathis 1996) to 
construct dust models that use up minimal amounts of carbon. The lower C/H ratio in the local 
ISM is marginally consistent with the lowest value of ~ 2.7 10 -4 present in PAHs, graphite, and 
in the form of C + that is implied by our observations. 



5.3. Comparison to the General Interstellar Extinction 

An additional constraint on the dust model is provided by the interstellar extinction. It is 
reasonable to expect that the extinction provided by the dust that fits the average high-latitude 
ISM emission spectrum will also provide a reasonable fit to the average interstellar extinction 
curve. Figure 4a compares the extinction from model A with = 100 to the average interstellar 
extinction curve of Mathis (1990). For sake of the comparison, we have multiplied the observed 
extinction curve, which is normalized to the total H column density, by the Nh/Nhi ratio, which 
is equal to ~1.2 (Savage & Mathis 1979). The figure shows that the predicted extinction is similar 
to the average interstellar extinction for the H I gas. At wavelengths ^ 20 ^m the value of the 
extinction is uncertain by at least a factor of two (Mathis 1990), so that the model extinction is 
essentially consistent with the observations. 

At UV wavelengths differences between the predicted and the average observed interstellar 
extinction curve may reflect fundamental uncertainties in dust models, or may suggest a real 
excess in the high latitude extinction compared to the average Galactic value which is biased 
towards lower latitudes. The extinction predicted by the model is higher by a factor of ~ 1.3 
at A = 0.125 /im, increasing to a factor of 2 at 0.1 /xm (Figure 4b). The contribution of PAHs 
to the interstellar extinction is quite uncertain flVerstraete fc Leger 1992 , loblin, Leger, fc 



Martin 1992). The UV properties of the PAHs derived by Desert et al. were tailored to fit the 
observed interstellar extinction curve given in Savage & Mathis (1979). As originally suggested by 
Greenberg & Chlewicki (1983), the carriers of the 2200 A extinction feature in the Desert et al. 
model have an otherwise flat extinction curve, and represent a distinctly different population of 
particles from those responsible for the rise in the FUV extinction. A similar approach was taken 
by Siebenmorgen & Kriigel (1992), who adopted a small graphite grain component in their dust 
model, but without the characteristic rise in extinction at UV wavelengths (Fig lb in their paper), 
attributing it entirely to PAHs. 
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The excess of high latitude UV extinction predicted in this paper is therefore directly the 
result of the enhanced abundance of PAHs derived in our model (a factor of two larger than that 
derived by Desert et al.) and the adopted PAH UV properties. In our model, the 2200 A feature 
is provided by small graphite particles, which also contribute to the rising extinction in the FUV 
(see also Draine & Lee 1984, Figure 7). Clearly the UV properties of the PAHs should be modified 
if graphite grains are allowed to contribute to the UV extinction as well. A reduction in the rise 
of the UV absorptivity of the PAHs will reduce the excess UV extinction exhibited in Figure 4b. 

There is some observational evidence suggesting that the UV properties of the PAHs may 
need to be modified. Boulanger, Prevot, & Gry (1994) have searched for a correlation between 
the strength of the 12 and 25 /an IRAS emission and UV extinction from select lines of sight 
through the Chamaeleon molecular cloud complex. They find a correlation between the intensity 
of the 12 fim emission, which is a measure of the PAH abundance, and the strength of the 2200 A 
extinction feature, suggesting a common carrier. Furthermore, they found no correlation between 
the intensity of the 12 /j,m emission and the rise in the extinction towards UV wavelengths. These 
observations suggest that the UV properties of the PAH are quite different from those adopted 
here. 



6. THE FAR-IR SPECTRUM AND THE PRESENCE OF "COLD DUST" 

Dust model parameters were solely derived from fitting the DIRBE 3.5 — 240 /tm data, 
and Figure 3 shows the model spectrum extended to FIRAS wavelengths. The figure shows 
that even without any attempts to fit the FIRAS data, the model produces a good agreement 
with the 240 — 500 /tm spectrum. Figure 5 shows the fractional residual emission, defined 
as [F (FIRAS) — F (model)] /F (model), as a function of wavelength. The residual exhibits a 
systematic trend in the 100 — 500 /im wavelength region, consistent with an underestimate in the 
calculated dust temperatures. To demonstrate the sensitivity of the far-IR spectrum to changes in 
dust temperatures, we plotted in the figure the fractional residual expected if we had attempted 
to fit a spectrum with an intrinsic temperature of 17.5 K with a model at 16 K. The figure shows 
that such a difference could account for the residual emission in the 240 — 500 /tm wavelength 
regime. Such a temperature increase would, for example, require the intensity of the LISRF to be 
increased by a factor of 1.5, or alternatively, require the ratio of the UV-optical to IR emissivity 
of the dust to be increased by the same factor. 

Persistent in the observations is, however, an excess of emission at wavelengths ^ 500 /im 
above the calculated model spectrum. In each of the FIRAS channels the excess is only aid 
detection, but considering the similarity of the trend in all the 500 — 1000 /im channels, the 
excess gains significance. It was previously detected in the FIRAS data by Wright et al. (1991), 
Reach et al. (1996), and may also be present in the residuals between the fit of a dust spectrum, 
characterized by a single-temperature of T=17.5 K and a u 2 emissivity law, to an Hl-correlated 
DIRBE and FIRAS Galactic emission spectrum (see Boulanger et al. 1996, Figure 3). The nature 
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of this excess emission component is still uncertain. We emphasize that this emission component 
should not be confused with the much stronger one found at high latitude by Reach et al. (1995), 
which is claimed by Puget et al. (1996) to be of extragalactic origin. 

Assuming that the excess is not an instrumental effect, it should be of a Galactic origin, since 
it was derived by correlating the FIRAS data with a template of Galactic 100 emission. As 
a result, any extragalactic component such as the one reported by Puget et al. (1996) would 
not be present in this data. Several explanations have been suggested for a Galactic origin of 
this component QWright et al. 1991| , |Wright 1993] , Reach et al. 1995|) : (1) the excess could be 



due to a far-IR feature in the dust emissivity; or (2) it could be due to a Galactic population of 
cold (T ~ 6 K) dust particles. Their low temperature is not likely the result of shielding from 
the interstellar radiation field, requiring enhanced dust cooling mechanisms such as expected for 
fractal dust grain, or needles. Recent laboratory results on the millimeter absorptivity of silicate 
grain materials ( Agladze et al. 1996| ) show that the opacities of these grains is significantly larger 



than obtained by extrapolating the Draine-Lee values with a A emissivity law. 



7. SUMMARY 

We presented the COBE 3.5 - 1000 //m spectrum of the diffuse ISM at latitudes \b\ > 45° 
(Fig. 1) and showed that the near-IR intensity ratios cannot be produced by any grain size 
distribution of bare silicate or graphite grains (Fig. 2). In particular, we find that the 3.5-to-4.9 
lim. intensity ratio is consistent with PAHs being the dominant emission component at these 
wavelengths. We used the DIRBE data to derive an interstellar dust model that consists of PAHs 
and bare silicate and graphite grains, stochastically heated by the local interstellar radiation field 
(LISRF). In addition, we used the [C II] 158 ^m line observed by the FIRAS to derive the C + /H 
abundance in the diffuse high latitude ISM. The results of this paper can be briefly summarized 
as follows: 

1) The dust model provides a very good fit to the 3.5 — 1000 /im diffuse ISM spectrum 
(see Figure 3). Dust parameters are summarized in Table 3. The model requires essentially all 
the available Mg, Si, and Fe, and ~ 15 % of the available O to be locked up in silicate grains. 
About 20% of the total carbon is locked up in PAHs (C/H = (7.3 ± 2.2) 10~ 5 ) and about 60-70% 
is locked up in graphite (C/H ~ (2.5 ± 0.8) 10 ). The abundance of PAHs in our model is higher 
by a factor of ~2 compared to that derived in previous models of low latitude cirrus. 

2) The energy radiated by the PAHs per unit H atom is 1.6 10~ 31 W/H atom, and the 
energy absorbed from the LISRF per unit carbon atoms in PAHs is 2.6 10~ W/C atom. The 
ratio between the two quantities provides an independent estimate of the C/H ratio locked up 
in PAHs. The resulting value, (C / H)pah = 6.1 10 -5 , in very good agreement with the PAH 
abundance required to account for the IR emission. 
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3) The detection of the [C II] 158 /im line by the FIRAS provides important constraints 

on the gas phase abundance in the cold neutral medium (CNM). We estimate the C + /H ratio to 
be between ~ (0.5 — 1) 10~ 4 (see Appenix A. 3), constituting ~ 10-20% of the cosmic abundance 
of carbon. 

4) The total amount of carbon in PAHs, graphite grains, and in the gas is C/H ~ (4.0 ± 1.2) 10~ 4 . 
This value is consistent with the cosmic carbon abundance, but barely consistent with the 
interstellar medium value of C/H = (2.25 ± 0.5) 10 -4 . 

5) The model provides a good fit to the FIRAS spectrum in the 240 to 1000 /im wavelength regime. 
Detailed comparison of the residual spectrum (Figure 5) shows a systematic trend consistent 
with an underestimate of the interstellar dust temperature. The figure also shows the presence of 
an excess of emission at A ^ 500 /im above that predicted by the model. This excess emission 
correlates with the template of Galactic 100 /im emission, and is therefore of Galactic origin. The 
nature of this excess emission is still unresolved. 

6) Our model provides generally a good fit to the average interstellar extinction curve (see 
Figure 4a). However, it produces an excess of extinction at FUV wavelengths (Fig. 4b). This 
excess is mainly due to the higher abundance of PAHs in the model, and the fact that we used 
Draine-Lee graphite particles with rising extinction in the FUV in combination with PAHs which 
have similar optical properties at FUV wavelengths. A fraction of the excess extinction may be 
real, reflecting a trend of increasing UV extinction with latitudes. The trend may be the result of 
an increasing PAH abundance at high Galactic latitudes. 

J. Ballester acknowledges the support of the NASA JOVE Program. E. D. acknowledges 
the hospitality of the Institut d'Astrophysique Spatiale in Orsay where the final revision to the 
manuscript was made. We also acknowledge helpful conversations with Xavier Desert, Francois 
Boulanger, Bruce Draine, Derek Massa, John Mathis, and Jean-Loup Puget. 

A. CHARACTERIZATION OF THE DUST AND PAH MOLECULES 

A.l. Heat Capacities 

The heat capacity plays a special role in determining the spectrum of stochastically heated 
particles. Small values of C(T) will accelerate the cooling around T, reducing the flux emitted 
at the appropriate wavelengths. In contrast, large heat capacities, especially at very low 
temperatures, may result in an excess of long-wavelength emission. 

At sufficiently high temperatures, the heat capacity of the various PAHs can be calculated 
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using the group additive method described by Stein, Golden, and Benson (1977; hereafter SGB). 
In this method, a PAH is decomposed into four carbon groups characterized by the composition of 
their nearest neighbors. For simplicity, we designate these atomic groups as groups A — D. Group 
A consists of the CH group in benzene (CqHq); group B of unsubstituted carbon in naphthalene 
(CioHq); group C is found in phenanthrene (Cu-ffio), and group D consists of interior carbon, 
identical to that encountered in a monolayer of graphite . For example, chrysene (CisH^ consists 
of 12 members from group A, 2 from group B, and 4 from group C (it has no interior carbon 
atoms; see, for example, Fig. 22 of Allamandola, Tielens, & Barker (1989). The number of 
members from each group in a given PAH are not independent, and one can derive several simple 
relationships between them and N c = Na + Nb + Nq + Nd, the total number of carbons in a 
PAH. For a fully hydrogenated PAH, the number of hydrogen atoms is given by Nh = y/6N c 
(see also Omont 1986[ ), and Na = Nh from the group definition. Empirically, we find that 



Nb + Nc = Na — 6. To these four groups we added a group A* to represent a dehydrogenated A 
group, so that Na* = (1 — Jh)Na, where fn is the number of filled hydrogen sites (/# = for a 
completely dehydrogenated PAH). 

The heat capacity of a given PAH at temperature T, Cpah(T), is then given in the group 
additivity technique by the sum: 

C PAH (T) = f H N A C A {T) + (1 - f H )N A C A * (T) + N B C B {T) 

+N C C C {T)+N D C D (T) 1 ) 

where the N's are the number of group members in the PAH. Following SGB, we assumed that Cq 
is equal to Cb, and substitution of the values of the A^'s in the equation for Cpah(T), gives for a 
completely hydrogenated PAH (/# = 1) 

Cpah{T) = N c { ^ [C A + C b - 2C d ] + ^-[C D - C B ] + C D } (A2) 

The equation shows that for sufficiently large values of N c , the heat capacity becomes equal 
to the graphitic limit. The heat capacity of the various groups in the 300 to 3000 K temperature 
range is given by Stein (1978) and Stein & Fahr (1985). For Ca* we used the heat capacity of 
phenyl radical [designated as Cb — (•) by Stein and Fahr (1985)]. 

The resulting values for Ca, Cb, and Cb can be written as polynomials in temperature: 

6 

C(J mole- 1 R- 1 ) = ^ a n T n , (A3) 

n=0 

where for C A : 

{a _ 6 } = {-1.23, 4.9 10" 2 , 2.07 10" 5 , -6.93 10~ 8 , 4.85 10 -11 , -1.44 10" 14 , 1.57 10~ 18 } 
for C B : 

{a - 6 } = {-7.02, 48.87 10" 2 , -1.11 10~ 4 , 7.98 10~ 8 , -3.33 lO" 11 , 7.42 10~ 15 , -6.79 10~ 19 } 
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and for Ca* ' 

{a _ 6 } = {4.01, 1.62 10~ 2 , 4.89 1(T 5 , -1.20 10~ 7 , 1.07 lO" 10 , -4.41 10" 14 , 6.91 lO" 18 } 

Since the group additive method breaks down at PAH temperatures below ~ 300 K, the 
heat capacity at these temperatures was fitted with that of graphite, scaled to smoothly merge 
with that calculated by the group additive method. The resulting value for the heat capacity in 
erg K , valid for T < 2000 K, is given by: 

log w [CpAH(T)/N c ] = -21.26 + 3.1688 log w T - 0.401894 (log w T) 2 (A4) 

In the limit of small carbon atoms, the PAH heat capacity reproduces tabulated values of 
naphtalene (C\oH 8 ), for which C(T) = 8Ca(T) + 2Cb(T), and anthracene (Cu-ffio), for which 
C(T) = 8Ca(T) + 2Cb(T). For large enough PAHs, the group additive method matches the 
heat capacity of graphite. The value of Cpah(T) in eq. (A4) is essentially identical to the one 
presented by Puget & Leger (1989). 

For graphite, the heat capacity in the 10 — 2000 K temperature is given in an easily integrable 
form by ( [Markela, Volga, fc Buchnev 1973| ): 

5 5 

Cgr f (T) = 2.2 x 10 7 a n T n / £ b n T n [ erg K^cmT 3 ] (A5) 

n=3 n=0 

where the coefficients are given by: {a 3 ,a 4 ,a 5 } = {0.10273, 4.4354 10" 2 , 2.2124 10~ 4 }, and 
{bo,bi,b 2 ,b 3 ,b4,b 5 } = {1.0003 10" 12 , 3.6909 10 4 , 1129.71, 30.4420, 1.2888 10" 2 , 1.0000 10~ 4 }. 

Silicate heat capacities were adopted from the values given by Draine & Anderson (1985). 
Silicate and graphite heat capacities were corrected by a factor of (1 — 2/N), where TV" is the 
number of atoms in the grain, to take into account that the heat from photon absorption events 
gets only distributed in the vibrational modes of the solid ( Guhathakurta fc Draine 1989| ). No 
such correction was needed for the PAHs since in the small particle limit their heat capacities 
reproduce that of observed molecules. 



A. 2. Cross Sections 

The cross section for graphite and silicate dust particles is given by: 

a(a,\) = ira 2 Q(a,\) (A6) 

where Q(a, A) is the dust absorption efficiency. Values of Q(a, A) were calculated from Mie theory 
using the dust optical constants from Draine & Lee (1984). 

The PAH cross sections used in this paper were adopted from Desert, Boulanger, & Puget 
(1990). In their model, the PAH cross section (JpahW is represented by a sum of UV-visual, IR 
continuum, and IR line terms, as follows: 

<7pah(A) = auv-visW + viRc(X) + o-/ra(A) (A7) 
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where 

auv-visW = W- 18 N c \ Pl f v (x)+ P2 f u (x)]C(£- c ) cm 2 
°iRcW = 3 -^^N c exp(-^) cm 2 (A8) 



(A- A,) 



2 



N c is the number of C-atoms in the PAH, x = l/A(/xm), p\ = 4.0 and p 2 = 1.1 are numerical 
parameters, x c = — 12,5 u , apAff(A) = 10-y/iV c /120 is the radius of a PAH, and A m = 10 /xm. 

The functions / u , /„, and C(y) are given by: 

f u {x) = (x-5.9) 2 (0.1a; + 0.41) for x > 5.9 /im" 1 

= otherwise (A9) 

/^(x) = 1.0 for x > x\ = 4 /im _1 

= x 2 (3x; — 2x)/xf for x <X[ 



C(y = x/x c ) = vr" 1 arctan(10 3 (y - l) 3 /y) + 0.5 (A10) 

The last expression in eq. (A8) assumes that the PAH spectral features are Gaussian in 
shape. The cross section is represented as a sum over the PAH feature j, where tjj is the value of 
the cross section at the line center Xj, {o^/lO^ 21 cm 2 }={35N H , 4.1N C , 2.9N C , 3.0N H , 47N H }, 
for, respectively the 3.3, 6.2, 7.7, 8.6, and 11.3 /im. bands (Leger, d'Hendecourt, & Defourneau 
1989; hereafter LHD). The AAj's are the effective width of the j features, satisfying the condition 
that the integral of the cross section over the line equals Aj = Xj(/j,m) 2 x /j/1.13 x 10 20 , where 
fj is the dimensionless oscillator strength (LHD). The linewidths in our representations are then 
given by: A.Xj(fj,m) = 3.53 x 10~ 17 fjXj(/j,m) 2 /<jj(cm 2 ). The resulting values for AXj(fim) are 
{0.0165, 0.070, 0.274, 0.157, 0.115} at {A j (^m)}= {3.3, 6.2, 7.7, 8.6, 11.3}, respectively. Note 
that our values for AXj(fim) are smaller by a factor of \^2tt from those quoted by LHD because 
our choice of the functional representation of the features. The actual functional representation 
is irrelevant for broad band photometry (such as the DIRBE data) as long as the condition that 
Aj = J aj(X)dX across a feature is satisfied. 



A. 3. The Carbon Gas Phase Abundance in the Neutral Medium 

The average high-latitude ISM spectrum calculated in the model represents the IR emission 
integrated over all gas phases along the line of sight (see eq. 1), including in addition to the H I, 
any high latitude molecular or ionized gas that correlates with the 100 fim emission. However, we 
do not expect emission from molecular gas to be a significant contributer to the IR spectrum from 
the average ISM gas. The high-latitude CO survey of Magnani, Blitz, & Mundy (1985) found that 
the surface filling factor of the high-latitude molecular gas is low, with estimates ranging from ~ 6 
to 10% (see review by [Magnani 1994 ). Rarely, molecular gas may be the dominant phase in select 
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cirrus clouds, with MBM 53-55 an extreme example (see Dwek, Arendt, Fixsen, D. J., & Reach 



1996). 



On the other hand, the warm ionized medium (WIM; a low density ionized gas that gives rise 
to pulsar dispersion measures, and diffuse Ha emission; e.g. Reynolds 1990), can be a significant 
contributor to the observed IR emission. Pulsar dispersion measures show that at high-latitudes, 
the 77 + /H I column densities are ~ 0.2 — 0.6 (e.g. Reynolds 1991). Assuming that the dust 
in the WIM is exposed to a radiation field similar to the LISRF and have a similar dust-to-gas 
mass ratio as in the general ISM, we get that at most one third of the observed IR intensity can 
originate from the ionized gas. 

We therefore assume that most of the IR emission originates from the neutral medium, and 
that almost all of the gaseous carbon in the average high-latitude ISM is in the form of C + , which 
dominates the contribution of heavy ions to the electron density. With an ionization potential of 
11.26 eV it is easily ionized by the diffuse interstellar radiation field, which is not hard enough to 
eject a second electron. C + has a ground-state fine-structure line at 158 /mi, which is collisionally 
excited by electrons, H I and H2. The amount of gaseous C can therefore be estimated from 
7(158 /im) observations. The 158 /im intensity can be calculated from the FIRAS observations 
(see Figure 1). Subtracting the dust continuum intensity calculated by the model we derive a 
value of 47r/(158 /j,m)/H = 1.45 x 10~ 33 W/H atom, a value lower by about a factor of ~ 2, 
compared to the |6| > 15° line intensity given by Bennett et al. (1994). 

Most of the contribution to the observed 1(158 /mi) intensity comes from the neutral gas 



( Bennett et al. 1994 ), instead of the WIM. The argument presented by Bennett et al. was simple: 
any [C II] 158 /mi emission from from the WIM will be accompanied by Ha emission. The relation 
between the two quantities is given by (Reynolds 1992): 

f + / 17+ 

7(158 /xm) ~ 1.45 x l£-57(_A_) I(Ha) (All) 

where 

I(H a ) ~ 8.7 x 10" n r 4 " - 92 EM W m~ 2 sr' 1 (A12) 

and where T 4 = T/10 4 K, and EM = J n\ds is the emission measure in units of cm 6 pc. Values 
of EM along several high latitude lines of sight are typically ~ 1 — 2 cm -6 pc, giving upper limits 
for the 7(158 /im) intensity of ~ (2 — 3) x 10~ 10 W mT 2 sr^ 1 , assuming that all the carbon is in 
the gas phase, and singly ionized. This intensity is significantly below the FIRAS observed value 
of ~ 2 x 10~ 9 W m~ 2 sr^ 1 flWright et al. 1991 , Bennett et al. 19~94]) . A similar conclusion was 



reached by Bock et al. (1993), who derived a value of 7(158 /im) (9 ± 3.5) x 10 U W m 2 sr 1 
for the line intensity from the WIM by extrapolating the 7(158 /an) versus Nh i correlation to 
N H i = 0. 

Since most of the 7(158 /im) emission originates from the neutral medium, we can can use 
the observed line intensity to derive an estimate of the [C + ]/[77 7] ratio in this medium, if the 
physical conditions in the medium (temperature, density, ionization state) are known. Designating 
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the upper and lower energy levels of the 158 /xm transition as levels 2 and 1, respectively, the line 
intensity (in erg s~ 1 cm~ 2 sr~ 1 ) is simply the integral of the line cooling rate along the line of sight 
(e.g. Osterbrock 1989| ), and can be written as: 



= ^ Nc+ ( 9 -l)e- h ^ kT E S'n-, (A13) 

where hv = 1.26 W~ 14 erg is the energy of the transition, N c + is the column density of C + 
atoms along the line of sight, A21 = 2.36 10~ 6 s" 1 is the Einstein A-coefficient for the 2 — > 1 
transition, 52 = 4 and g\ = 2 are, respectively the statistical weights of levels 2 and 1, rij is 
the number density of the colliding species, and n c j = ^42i/<?2i is the critical density, where 
Q21 = 8.629 10~6 Q,j / g2T l / 2 crm? s _1 is the collisional de-excitation rate, and f2j is the dimensionless 
collision strength with species "j". For electronic collisions, Q ~ 2.8 in the ~ 5000 to 20,000 
K temperature interval, and O ~ 1.8 for temperatures below ~ 1000 K (Hayes & Nussbaumer 
1984). For collisions with H, Q, = 0.00292 (Launay & Roueff 1977). Normalizing the observed 
line intensity to the observed H I-column density along the line of sight, eq (13) can be used to 
calculate the C + -to-H I abundance ratio, giving: 

[C+] _ (91/92) AttI(u) m 

[HI] A 21 hu [ N HI U Y n J 1 j 

where we have made the approximation that rij <C n c j. 

Most of the line emission originates from the cold neutral medium (CNM) which has typical 
parameters {ra#,T} = {20 — 80 cm -3 , 40 — 100 K}, with pressures P/k ps 2500 cm~ 3 K (e.g. 
Heiles 1994). This medium is kept partially ionized by high-energy cosmic rays (CR). The electron 
density, n e maintained by a CR ionization rate Qcr (hi s ) is given by QSpitzer 1978j) : 

n, = U i [1 + (1 + %^ 2 ] (A15) 

where m 6.2 10~ n T" 1 / 2 cm 3 s 1 is the partial H— recombination coefficient to the n>l levels 
of hydrogen, and rtj is the total number density of photo-ionized heavy elements. 

Equations (A14) and (A15) can be solved for [C + ]/[i7 I] as a function of the parameters of 
the CNM, the CR ionization rate CcRi an d the observed I(158/im) /Nh i ratio. Estimates 
of Ccr range from ~ 10~ 16 — 10~ 18 For the C + cooling rate per H-atom of 

47r I (158 nrn)/Njii = 1-45 10 -33 W/H atom found earlier, and an adopted pressure of 
P ps 2000 - 3000 cm- 3 K, we get a [C+]/[HI] ratio of: 

jgj « (0.5 - 1) x 10- 4 (A16) 

For a mean atomic weight of \i = 1.42, the resulting C + -to-gas mass ratio, Z c + ~ (4 — 9) x 10~ 4 . 
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Table 1. Dust Spectra/color temperatures in the Average ISM and Selected Cirrus Clouds 3, 



1(A) //( 100 /iin) / T color (K) 



A(/xm) 


average ISM b 


Cloud l c 


Cloud 2 d 


3.3 


0.00168 


0.00161 


0.00124 


4.9 


0.00252 /550 


0.00074 /1500 


0.00186 /500 


12 


0.0464 /230 


0.0343 / 200 


0.0600 /210 


25 


0.0389 /210 


0.0283 / 210 


0.0633 /190 


60 


0.160 / 70 


0.114 / 70 


0.178 / 75 


100 


1.0000 / 23 


1.0000 / 44 


1.000 / 24 


140 


1.93 / 18 


2.53 / 16 


2.45 / 16 


240 


1.28 / 19 


2.51 / 16 


1.36 / 21 



a Fluxes are color-corrected and given in units of 0.7 MJy sr~ l /Nh /(10 20 cm 2 ). Numbers in 
parethesis in each row "j" are the color temperatures determined from the I(Xj-i) / 'l(Xj) flux ratios 
for a A~ 2 dust emissvity law 

b Fluxes represent Galactic average IR fluxes integrated over all gas phase components that 
spatially correlate with the observed H I emission (see §2 for more detail) 

c average spectrum of the cloud complex MBM 53, 54, and 55 (see §2 for more details) 
d average spectrum of cloud A (Low et al. 1984) 
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Table 2. Dust Model Parameters 


Parameter 


Value 






Model A 


Model B 


PAH parameters 






N cl 


20 


20 


N c2 


variable 


variable 




variable 


variable 


r/ 

^PAH 


variable 


variable 


graphite parameters 


model A 


model B a 


a min (fim) 


variable 


1.29 x W~ 4 Nll 3 


71 


-3.5 


3 7p + 2 


a b ([im) 


0.1 


variable 


72 


-3.5 


-3.5 




0.25 


0.25 


Z grf 


variable 


variable 


silicate parameters 






a min (pm) 


0.0050 


0.0050 


7 


-3.5 


-3.5 


a max (nm) 


0.25 


0.25 




variable 


variable 



a In model B, a m i n is calculated so that the number of carbon atoms in a spherical graphite grain 
of that radius is equal to N C 2, the number of carbon atoms in the largest PAH, and 71 is calculated 
so that dn(N c ) / dN c is continuous across a m i n , the PAH-graphite radius boundary 
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Table 3. Summary of Dust Model Parameters 3, 



Model Parameters 


Averaj 

Model A 


;e Diffuse ISM Dust 

Model B 


PA H 










1 v ci 


20 


20 


20 


20 




100 


171 


100 


163 


7 P 


-1.67 


-2.33 


- 1.78 


-2.67 


ZpAH 


0.00062 


0.00061 


0.00060 


0.00054 


\jr i UjJJI oboe 












0.0005 


0.0005 


0.0006 


0.0007 


7i 


-3.5 


-3.5 




-6.00 


a b ([im) 


0.10 


0.10 


0.0006 


0.0012 


72 


-3.5 


-3.5 


-3.5 


-3.5 




0.25 


0.25 


0.25 


0.25 


Z grf 


0.0021 


0.0020 


0.0022 


0.0023 












a min (fim) 


0.0050 


0.0050 


0.0050 


0.0050 


7 


-3.5 


-3.5 


-3.5 


-3.5 




0.25 


0.25 


0.25 


0.25 




0.0053 


0.0058 


0.0046 


0.0042 


SWid Phase Abundances 










^carbon 


0.0027 


0.0026 


0.0028 


0.0028 


Zdust 


0.0080 


0.0083 


0.0074 


0.0070 


Gas Phase Carbon Abundances 










z c+ 


0.0004-0.0009 


0.0004 


-0.0009 



a See Table 2 for the definitions of Models A and B 

b Dust-to-gas mass ratios are given in units of jf^j- For comparison, the cosmic abundances of 
refractory elements that can be locked up in dust are: Z car f) 0n (Q) = 0.0030; Z s h(Q) = 0.0043 
assuming a silicate composition of {MgSiFe}0^. Elemental abundances were taken from Anders 
& Grevesse (1989). 

c Calculated from the 7(158 /xm)//(100 jum) flux ratio. See Appendix A. 3 for more details 
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Fig. 1. — The various dust spectra observed by DIRBE and FIRAS (the latter is given only for the 
average ISM). The average spectrum of dust in the diffuse ISM normalized at 100 /um to a flux of 
0.7 MJy sr^ 1 /Nh /(10 20 cm -2 ). The other spectra have arbitrary normalizations, and are offset 
from the diffuse ISM spectrum for sake of clarity. 



Fig. 2. — A comparison of observed intensity ratios with calculated DIRBE intensity ratios for 
PAHs and individual grains of various sizes and compositions. Intensity ratios for the individual 
PAHs were calculated as a function of N c , the number of carbon atoms in the molecule, and for 
illustrative purposes converted to a radius using the relation: a(A)=0.913y/N c . 



Fig. 3. — Fit of our interstellar dust model to the average dust spectrum in the diffuse ISM for 
model A with N& = 100. Details of the grain size distribution and relative abundances of the 
various dust components can be found in the text (§4, and Table 3). 



Fig. 4. — The average extinction for the diffuse interstellar medium derived from model A with 
N C 2 = 100 (see Table 3), is compared with the observed average interstellar extinction curve 
presented by Mathis (1990) renormalized to H I column densities. Detailed of the comparison are 
presented in the text. 



Fig. 5. — Residuals of between the FIRAS average ISM spectrum and model calculations. The 
residuals show a trend in the 100 — 500 /im region consistent with an underestimate in the 
temperature of the emitting dust. The smooth curve depicts the residuals expected with such 
an underestimate (see text for details). However, even with such a correction, the figure shows a 
persistent excess of emission at A > 500 (im. 



